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Abstract — n this paper we show how appropriate 
superpositions of Bessel beams can be successfully used to 
obtain arbitrary longitudinal intensity patterns of nondiffracting 
ultrasonic wavefields with very high transverse localization. 
More precisely, the method here described allows generating 
longitudinal acoustic pressure fields, whose longitudinal intensity 
patterns can assume, in principle, any desired shape within a 
freely chosen interval < z < L of the propagation axis, and 
that can be endowed in particular with a static envelope (within 
which only the carrier wave propagates). 

Indeed, it is here demonstrated by computer evaluations that 
these very special beams of non-attenuated ultrasonic field 
can be generated in water-like media by means of annular 
transducers. Such fields "at rest" have been called by us Acoustic 
Frozen Waves (FW). 

The paper presents various cases of FWs in water, and 
investigates their aperture characteristics, such as minimum 
required size and ring dimensioning, as well as the influence 
they have on the proper generation of the desired FW patterns. 
The FWs are particular localized solutions to the wave equation 
that can be used in many applications, like new kinds of devices, 
such as, e.g., acoustic tweezers or scalpels, and especially 
various ultrasound medical apparatus.n this paper we show how 
appropriate superpositions of Bessel beams can be successfully 
used to obtain arbitrary longitudinal intensity patterns of 
nondiffracting ultrasonic wavefields with very high transverse 
localization. 

More precisely, the method here described aUows generating 
longitudinal acoustic pressure fields, whose longitudinal intensity 
patterns can assume, in principle, any desired shape within a 
freely chosen interval < z < L of the propagation axis, and 
that can be endowed in particular with a static envelope (within 
which only the carrier wave propagates). 

Indeed, it is here demonstrated by computer evaluations that 
these very special beams of non-attenuated ultrasonic field 
can be generated in water-like media by means of annular 
transducers. Such fields "at rest" have been called by us Acoustic 
Frozen Waves (FW). 

The paper presents various cases of FWs in water, and 
investigates their aperture characteristics, such as minimum 
required size and ring dimensioning, as well as the influence 
they have on the proper generation of the desired FW patterns. 
The FWs are particular localized solutions to the wave equation 
that can be used in many applications, like new kinds of devices, 
such as, e.g., acoustic tweezers or scalpels, and especially various 
ultrasound medical apparatus.I 
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I. A Preliminary Introduction 

THE theory of Localized Waves (LW), also called Non- 
Diffracting Waves, has been developed, experimentally 
verified and generalized over the years. The application sectors 
may include diverse areas such as Optics, Acoustics and 
Geophysics. One of their striking properties is that their peak- 
velocity can assume any value between zero and infinity Q. 
The most important characteristic of the LWs, however, is that 
they are solutions to the wave equations capable of resisting 
the effects of diffraction, at least up to a certain, long field- 
depth L. They exist both as localized beams and as localized 
pulses. The LWs that have been more intensely investigated 
are the so-called "superluminal" ones|2|, they being the easiest 
to be mathematically constructed in closed form, and therefore 
experimentally generated. Let us recall that the first LWs of 
this kind have been mathematically and experimentally created 
by J.-y. Lu et al. in Acoustics|3|, |4| (when they are actually 
supersonic, and not superluminal), and used by J.-y. Lu et al. 
for the construction of a high-definition ultrasound scanner, 
directly yielding a 3D image|5|. 

More recently, also for the subluminal LWs various analytic 
exact solutions started to be found|6|. Among the subluminal 
localized waves, the most interesting appear to be those 
corresponding to zero peak- velocity, that is, those with a static 
envelope (within which only the carrier wave propagates): 
Such LWs "at rest" have been called Frozen Waves (FW) by 
us|7l, [(Sl, f9|. The FWs have been actually produced for the 
first time, quite recently, in the sector of Optics ifTOl . 

The above information confirms that LW pulses and beams 
can be used also in ultrasound imaging; but this is not as 
evident for the FWs because of the static nature of their 
envelopes. What FWs make possible, by contrast, is the 
localization of arbitrary spots of energy within a selected space 
interval < z < L. 

Let us recall here some previous work done in connection 
with ultrasonic non-diffracting fields. This list, obviously, is 
not at all exhaustive, but is merely illustrative of previous 
literature related with pulses Ol, IH, iHTI . |[T2L with single 
Bessel beamsfTSl, fT4l, fTS^, with characteristics of the field 
emitted by flat annular arravs lT6l , ifTVl , and with scattering 
produced by spherical objects| lH, |[T9ll , ||2Q| . Many references 
therein could also be considered. Afterwards, a lot of work 
has been produced by investigating all the possible superpo- 
sitions of Bessel beams obtained via integrations over their 
axicon angles (that is, their speed) and/or frequencies and/or 
wavenumbers and/or phases, etc.: See, e.g., Q, (21 and the 
references quoted therein. 
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However, only few publications 1211 . ll22ll have addressed 
the superposition of Bessel beams with the same frequency 
but with different longitudinal wavenumbers: In particular for 
obtaining Frozen Waves (71, ISl, 191. 

Purpose of this paper is to contribute to the last topic, by 
the application of a general procedure previously developed 
by us for the generation of FWs (mainly in Optics), and then 
simulate the production of ultrasonic FW fields in water. 
Namely, we shall use the methodology in Refs.|7l, |8l, (91, 
that allows controlling the longitudinal intensity shape of the 
resulting fields, confined, as we were saying, inside a pre- 
selected space interval < z < L; where z is the propagation 
axis, while L can be much larger than the wavelength A of 
the adopted ultrasonic monochromatic excitation. In practice, 
we shall perform appropriate superpositions of zero-order 
(m = 0) Bessel beams. The generated fields, apart from 
possessing a high transverse localization, are endowed with 
the important characteristic that inside the chosen interval 
they can assume any desired shape: For instance, one or 
more high-intensity peaks, with distances between them much 
larger than A. 

Before going on, let us add a couple of preliminary com- 
ments about the FWs. The first is that it would be of course 
possible to use higher-order (m ^ 1) Bessel beams. In this 
case, however, it would be practically necessary a subdivision 
of the radiator ring elements into small arc segments (i.e., a 
segmented array fT^ ). because of the azimuthal phase depen- 
dence e^^^ of those fields on the beam axis z. It appears 
to be more convenient tackling with such complexities in 
the aperture design only after that an investigation of the 
cylindrically symmetric zero-order Bessel beam superpositions 
has been exploited. Anyway, we shall briefly come back, 
in Sec.2, to the higher-order Bessed beams question. The 
second comment is related to the flux of energy within a FW, 
especially in the realistic case of a finite aperture. Let us go 
back for a moment to a (truncated) Bessel beam, by recalling 
first of all that its good properties are due to the fact that, 
even in presence of diffraction, the "intensity rings" that are 
known to constitute its transverse structure (and whose values 
decrease when the spatial transverse coordinate p increases) 
go on reconstructing the beam itself all along a (large) depth 
of field. Namely, given a Bessel beam and a gaussian beam — 
both with the same energy E, the same spot Apo and produced 
by apertures with the same redius R in the plane z = — 
, the percentage of energy E contained in the central peak 
region < p < Apo is smaller for a Bessel rather than 
for a gaussian beam: It is just such a different distribution of 
energy on the trasnsverse plane that causes the reconstruction 
of the Bessel beam central peak even at large distances from 
the source, and even after an obstacle with sizes smaller than 
the aperture'^. Such a property is possessed also by the other 
localized waves. In other words, a certain energy must be 
contained in, and carried by, the side-lobes! Of course, such 
energy can be reduced towards its minimum necessary value 
by suitable techniques: See, e.g., Refs.|24|, |25|. 

Also the FWs we are dealing with in this paper are non- 
diffracting beams: and we investigate their generation by a 



limited aperture, so that their energy flux does remain finite. 
But they are just the energy contributions coming from the 
lateral rings that strengthen the FW field pattern, in the z 
direction, all along the FW field depth. 

The paper is organized as follows: After this Introduction, 
Section 2 briefly describes the methodology used for the 
generation of FWs. Then, in Section 3 the method for the 
computation of the ultrasonic fields is presented; followed in 
Section 4 by a discussion about the annular radiator require- 
ments. With respect to the last point, let us observe that, even 
if the precision requirements that one meets for the annular 
transducer sizes are tight, and appear as a real challenge to be 
overcome, nevertheless it is encouraging to know that various 
successes have been already obtained in similar situations, for 
example for monolithic|ll |, piezocomposite|23|, and PVDF 
transducers 1261 . More explicitly, let us mention that ultrasonic 
beams have been technologically synthesized, by superposing 
Bessel radiations, in interesting works like 1271 . 1281, 021, 

1291. (301, 1211, m, (13, (H. 

Finally, four simulated examples of ultrasonic FW fields 
in an ideal water-like medium assuming no attenuation are 
presented in Section 5; while some conclusions appear at the 
end of the paper. 

II. Introduction to the Frozen Waves (FW) 

In this Section, a method for the generation of FWs is pre- 
sented. For more details the reader is referred to reference|8|. 
The main aim here is constructing, inside the finite interval 
< 2; < L of the propagation axis {p — 0), a stationary 
intensity envelope with a desired shape, that we call \F{z)\^ . 
To such a purpose, we shall take advantage of the field 
localization features of the axis- symmetrical zero-order Bessel 
beams. Of course, we could start from higher order Bessel 
beams, and their relevant formalism. But, for simplicity, as 
already said, we prefere to start from the FWs given by the 
following finite superposition of zero-order Bessel function of 
the 1*^ kind, all with the same frequency but with different 
(and still unknown) longitudinal wavenumbers /^^ : 

N 

In Eq.([T]), quantities kpn are the transverse wave numbers of 
the Bessel beams, linked to the values of by the following 
relationship: 

^pn + /^n = 5 (2) 

where cjq = 27r/o, is the angular frequency, and c the plane 
wave phase velocity in the selected medium. By the way, 
let us recall that we leave the frequency fixed, since we are 
considering beams (and not pulses): We vary, afterward, the 
longitudinal wavenumber (and amplitude and phase) of each 
Bessel beam, so as to obtain the desired longitudinal shape of 
the resulting beam. 

It is important to notice that we restrict the values of the 
longitudinal wave numbers to the interval 
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in order to avoid evanescent waves and to ensure forward 
propagation only. 

Our goal is now finding out the values of the (3^^ and of 
the constant coefficients An in Eq[T] in order to reproduce 
approximately, inside the said interval < z < L, the desired 
longitudinal intensity pattern In other words, for p = 

we need to have: 



N 



AnC 



F{z) with 0< z <L . 



(4) 
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To obtain this, one possibility is to take f3n = ^f^, thus 
obtaining a truncated Fourier series that represents the desired 
pattern F{z). This choice, however, is not very appropriate 
because of two reasons: (i) it yields negative values for Pn 
when n < 0; which would imply backward propagations; (ii) 
we wish to have L ^ A, and in this case the main terms in 
the series could lead to very small values of the Pn, resulting 
in a very short depth of field, and strongly affecting, therefore, 
the generation of the desired envelopes far form the source. 
A possible way out is to take 

27T 



Pn 



(5) 



where the value of Q > can be freely selected. Then, on 
introducing Eq.([5]) into Eq.Q, the transverse wavenumbers 
can be expressed as 
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Now, inserting again Eq.(|5]i for /?„ into ([T]( and putting p = 0, 
one gets 
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Expressions ^ and ^ represent an approximation of the 
desired longitudinal pattern, since the superposition in ([T]) is 
necessarily truncated. However, we can control, and improve, 
the fidelity of the reconstruction by varying the total number 
of terms 2 + 1 in the finite superposition by a suitable choice 
of the parameters Q and/or L. The complex coefficients An in 
([8]) forward the final amplitudes and phases for each one of the 
Bessel beams in Eq.([T]); and, because we are adding together 
zero-order Bessel functions, we can expect a high degree of 
field concentration around p = 0. 

The methodology introduced here deals with control over 
the longitudinal intensity pattern. Obviously, we cannot get a 
total 3D control, owing to the fact that the field must obey 
the wave equation. However, we can get some control over 
the transverse spot size through the parameter Q. Actually, 
Eq.Q, which defines our FW, is a superposition of zero-order 
Bessel beams, and therefore we expect that the resulting field 
possesses an important transverse localization around p = 0. 
Each Bessel beam in superposition ([T]) is associated with a 
central spot with transverse width Apn ^ 2A/kpn\ so that. 



on the basis of the expected convergence of series ([T]), we can 
estimate the radius of the transverse spot of the resulting beam 
as being: 

2.4 

^ , = . . _ . (9) 
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Let us explicitly notice, at this point, that one can increase 
the control over the transverse intensity pattern of the resulting 
beam by having recourse to the already mentioned higher order 
Bessel beams. An ex.-p{—iujt)Ji^{kpnp) exp(z/3nz), with u > 
1, in the superposition ([T]), keeping the same values of the 
An given by Eq.([8|, and of the Pn given by Eq.Q. On doing 
this, the longitudinal intensity pattern can be shifted from p = 
to a cylindrical surface of radius p = p\ which can be 
approximately evaluated through the relation 
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This allows one to get interesting nondiffracting fields, mod- 
elled over cylindrical surfaces. Some details can be found in 
Ref.|8|. 

Anyway, relationship ^ itself is rather useful, because, 
once we have chosen the desired longitudinal intensity pattern 
we can also choose the size of the transverse spot 
(2Ap), and then use Eq.([9]) to evaluate the corresponding, 
needed value of Q. 

Such a method for creating FWs depends on the nature of 
the waves to be used. For instance, for an efficient generation 
in Optics, a recourse to ordinary lenses is required even when 
using the classical method with an array of so-called Dumin et 
al.'s circular slits [|35l , 1361 , 1371 . In several other cases, optical 
axicon lenses have been used |[38l , 1391 even if, to say the truth, 
axicon devices had been already proposed since long time|40| 
in Acoustics too. Of course, whenever possible, one may have 
recourse also to holographic elements, or suitable mirrors, etc. 
For instance, optical FWs were produced experimentally (for 
the first time) by using computer generated holograms and a 
spatial light modulator| 10|. 

In the case of ultrasound, we may adopt annular transducers 
composed of many narrow rings. Each one of them will have 
to be exited by a sinusoidal input having a particular amplitude 
and phase. Section 4 will discuss these points in more detail. 

III. Method for Calculating Ultrasonic Fields 

This Section summarizes the well known impulse response 
(IR) method, which is the mathematical basis used in this work 
for the computation of the ultrasonic frozen waves. Such a 
spatial impulse response method has recourse to the linear 
system theory to separate the temporal from the spatial char- 
acteristics of the acoustic field. Then, the aperture IR function 
can be derived from the Rayleigh-Sommerfeld formulation of 
diffraction, that is, from the Rayleigh integral lt41J . |i42J, Ii43i 
1.44.1 . via the following expression: 



L 



dS . 



(10) 



is 27r|ri-ro| 

Here the radiating aperture S is mounted on an infinite rigid 
baffle; quantity c being now the speed of sound, and t the 
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time. The spatial position of the field point is designated by 
ri, while tq denotes the location of the radiating aperture. 



The integral 10 is basically the statement of Huyghens 
principle, and evaluates the acoustic field, i.e., the field 
pressure relative to the unperturbed static pressure Pq, by 
adding the spherical wave contributions from all the small area 
elements that constitute the aperture. This process can also be 
reformulated by using the acoustic reciprocity principle, and 
then constructing h by finding out the part of the spherical 
waves that intersect the aperture. This implies that the IR 
function h depends both on the form of the radiating element 
and on its relative position w.r.t. (with respect to) the point 
where the acoustic pressure is calculated. The IR formulation 
assumes a flat or gently curved aperture (that is, the aperture 
dimensions are large compared to the wavelength), which 
radiates into a homogeneous medium with no attenuation]^ 
and operates in a linear regime. 

On using the spatial impulse response h{x,y , z^t), the final 
acoustic pressure can be written as 1421 . B31 : 

dv{t) 



p{x,y,z,t) = pr, 



dt 



* h{x,y,z,t) 



(11) 



where pm is the density of the medium, and the symbol star, *, 
denotes the convolution in time of the derivative of the velocity 
surface signal v, for a radiating element with aperture's IR 
function h. 

The above procedure was implemented into a matlab pro- 
gram, employing the DREAM toolbox B6ll , WT\ for the calcu- 
lation of the IR of the rings that compose the annular aperture. 
The steps followed by the program for computing the acoustic 
fields can be summarized as follows: 

1) Compute the theoretical FW desired pattern. 

2) Determine the transducer requirements. 
Sample the amplitude & phase FW profiles at 2; = 0. 
Assign sinusoidal Vk for ring Uk using sampled values. 
Sweep field points Xij;yij; Zij for transducer ring n/^. 



3) 
4) 
5) 
6) 
7) 
8) 



At point P{ij) calculate h and P(ij) = Pm^fc * ^• 



Accumulate pressure for ring Sl go back to step 4. 
Store and display results. 
The medium selected for the simulations was an ideal water- 
like lossless mediumj^ having a constant propagation speed 
of c = 1540 m/s [which is an average|48| among Cfat = 
1450 m/s, Cbiood = 1575 m/s, and Cmuscie = 1600 m/s]. 

IV. Aperture Dimensioning 

In this Section we discuss some dimensioning aspects of an 
annular aperture for creating the FW fields. They include: i) 
the selection of the transducer size; and ii) the number and 
dimensions of the rings. 

A. Transducer size 

The first element to be considered is the minimum aperture 
radius, (i^min), required for the generation of the desired 

^It is possible to include the attenuation of the medium by filtering the 
Fourier transform of /i by a frequency dependent attenuation function (421 . 
L45|. 

^The case including attenuation will be investigated in a forthcoming paper. 
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Fig. 1. Relationship between the aperture radius R and the summation 
limits entering Eq.|4}, in terms of the parameters L and wq = 27r/o. 
Continuous lines refer to frequency /o =2.5 MHz, while dashed lines refer 
to /o = 4 MHz. 



ultrasonic FW fields. The selection of this value, however, 
is not entirely free, because it depends on the chosen N, that 
determines the number i2N + 1) of Bessel beams added in 
Eq.Q. From Eqs.Q and ^ we can write: 

/3.ax = Q+^<^ (12) 



Q 



L 

27rN 



> . 



Afterwards, if we put /3max = ^ into Eqs.(12), and subtracts 
both sides of the equations, we get 



B ■ 

C 



47tN 



The equation for the maximum axicon-angle is thus 

, /J -^min kf)^^^ 

tan (7max — ^ 



(13) 



(14) 



Inserting the expression for kp from Eq.Q into Eq.( 14 ( and 
solving it for /^min, we get: 



(15) 
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Finally, is we equate expressions ([13] and [15]), we obtain for 
Rmin the expression 



-^min — La 



L } 



(16) 



which does allow the estimation of the minimum aperture- 
radius for the selected value of TV, with cjq and L as param- 
eters. 

If requisite (16) is not fulfilled, and R < Rmin, the resulting 
FW pattern may result distorted. This is because the higher 
is the precision desired for the FW, or the larger its desired 
maximum distance L, the higher will be the values needed for 
and the aperture radius R^. This is a logical conclusion, 
if we think the frozen waves as nothing but the realization of 
Huyghens principle for constructive/destructive interference. 
We can also get the above mentioned dependence by putting 
Pmin = in Eq.([T3]), and then solving it for N, so as to obtain 
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Fig. 2. Relationship between the aperture radius R and the number 
which determines the number of the Fourier terms entering Eq.JiJ, for a fixed 
distance L = 0.040 m but for various values of /q. 
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Fig. 3. Side view of the three simulated profiles, for a simple (ideal) FW 
consisting of two step-functions only [see case Case 1 of Sec. 5: Cf. Eq.(17) 
below], obtained by varying the working frequency. The parameter settings 
are the following: = 12; R = 35 mm; d = 0.3 mm; = 0.05 mm; 
Nr = 101; L = 30 mm, and /o = 2.5 MHz 



the maximum number of allowable terms A^max (this implies 
setting Q = in the Fourier superposition: 

N < TV^ax = ^ . (17) 
47rc 

However, care has to be taken when using using this limiting 
value, since usually it leads to unpractical transducer sizes 
(R > 50 mm), with too many rings (Nr > 100). 

Relationship ([16]) is at work in Figures [T] and [2] where one 
of the parameters (L or wq) is alternatively fixed, while the 
other is left variable]^ Both Figures clearly show how the value 
of R rapidly increases when N is increased. By contrast. 
Figure [T] depicts Eq.([T6]) for different values of L, using the 
two frequencies /o = 2.5 MHz, and /o = 4 MHz; and show 
how in both cases, up to a certain value of [e.g. N ^ 15 for 
the first case, continuous lines], the value of R almost does 
not change when the distance L is increased. 

The rapid increase of R with N can be partially mitigated 
by raising the working frequency /q , as shown in Fig [2] were a 
fixed distance L = 0.040 m is assumed. This option, however, 
has the effect of requiring smaller ring- widths (d^) for a good 
generation of the fields. 

In Fig|3]we can see how the frequency affects the simulated 
profiles, in the simple case of an ideal FW consisting of 
two step-functions only [a situation better exploited in Case 
1 of Sec. 5: Cf. Eq.(17) below]. The three intensity profiles 
simulated in Fig [3] have been obtained by varying the working 
frequency, always keeping the same aperture. Notice how, 
for frequencies larger than / ^ 2.5 MHz, any increase in the 
frequency distorts more and more the envelope of the desired 
FW. This effect is due to two causes: First, the dimensions 
of the rings are not changed, and , second, the patterns that 
result from the beam superpositions in Eq.([T]) are influenced 
by the approximately linear relationship existing between the 
ultrasound wavelength and the average distance among the 
sidelobe peaks of the superposed Bessel functions. At the 
end, the combination of these two causes produce a destructive 
interference effect on the generated FW. 

^We'd like to point out in these Figures that the estimated size of the 
spot-radius given by Eq.|9} changes when is varied. This is because Q 
disappears during the derivation of Eq.fT6). However, for the current setting, 
the changes are not significant and are in the range 0.3 < Ap < 0.5 mm. 
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Fig. 4. Comparison of of the simulated profiles, of the FW considered in 
Fig. 3, obtained now by using three different values for the transducer radius 
R. Settings: = 9; d = 0.3 mm; = 0.05 mm; L = 30 mm; 
/o = 2.5 MHz; A^^i = 115; Nr2 = 101, and Nrs = 86. 

The effect of changing the emitter radius R can be observed 
in Figure |4] which still refers to the FW considered in the 
previous graphic: That is, to an ideal FW consisting of two 
step-functions only [a case better exploited, let us repeat, in 
Figure [8] of Case 1 in Sec. 5: Cf. Eq.(17) below]. This time, 
the three simulated intensity profiles have been obtained using 
the different transducer radii Ri = 40 mm, R2 = 35 mm and 
Rs = 30 mm. Other parameters used for the simulation are: 
L = 30 mm, d = 0.3 mm, = 0.05 mm, N = 9 and 
/o = 2.5 MHz. Notice how, when the radius is reduced (e.g., 
R = 30 mm), the FW intensity profile gets more distorted 
with respect to the ideal envelope. 

B. Dimensions of the Rings 

The next step in the definition of the annular radiator is the 
dimensioning of the transducer rings. This means determining: 
(i) the width d, and (ii) the inter-ring spacing, or kerf, A^. 
On this respect, one meets a practical problem when the 
required distances are very small, for instance sub-millimetric 
dimensions. In this case, the technology employed in the 
annular transducer fabrication may strongly influence its per- 
formance. To illustrate this effect on the generated ultrasonic 
fields, in Figure [5] we show how small changes in d or 
A^ can lead to severe distortions in the profiles (and the 
same will be true, e.g., for Fig. 8 below). The corresponding 
parameters (d^; A^^.) chosen for this Figure are: First choice: 
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Fig. 5. Comparison of the simulated FW profiles when using different values 
ford; Ad. Settings: = 35 mm; = 9; L = 30 mm; /o = 2.5 MHz; 
Nri = 101; Nr2 = 88; A^^s = 88 and Nr4 = 117. 
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Fig. 6. Sampling of the ideal symmetric amplitude pattern, for the generation 
of the field in Fig[8] below. The dots indicate the sampled values used for the 
excitation of the rings. 
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Fig. 7. Sampling of the ideal symmetric phase pattern, for the generation 
of the field in the following FigjS] The dots indicate the sampled values used 
for the excitation of the transducer rings. 
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Fig. 8. The theoretical FW pattern chosen in the Case 1 of Sec. 5. Settings: 
N = 12; R = 35 mm; d = 0.3 mm; = 0.05 mm; A^^ = 101; 
L = 30 mm; and /o = 2.5 MHz. 



(0.30; 0.05) mm; Second choice: (0.35; 0.05) mm; Third 
choice: (0.30; 0.10) mm; and Fourth choice: (0.25; 0.05) 
mm. 

One can observe how, for dimensions larger than d ^ 
0.30 mm and A^^ ^ 0.05 mm, the FW envelopes gets heavily 
modified. This phenomenon is probably related to the above 
mentioned causes of destructive interference, as well as to a 
near-field effect of the radiator. One can attempt at partially 
overcome it by increasing the maximum distance of the 
FW fields. 

What discussed above leads to the following question: 
What are the dimensions of d and A^^ that will permit us 
to generate many different ultrasonic FW patterns with the 
same transducer? Stated in another way: What will be the 
ring dimensions (d; A^) that best comply with a predefined 
set of FWs to be generated? To address these questions, one 
has first to consider the sampling process involved during the 
generation of the frozen waves. This is a critical point in order 
to attain a proper excitation while minimizing the hardware 
cost. Indeed, each transducer ring will require in practice one 
electronic emission channel in order to operate]^ For instance, 
the number of rings (N^) and emission channels needed for the 
physical realization of the patterns shown in Figures [8}|T5] is 
in practice of A/"^ = 101; which means a considerable amount 
of electronic hardware for any application. In conclusion, 
besides the required electronics, the definition of the annular 

^It should be remembered that we are trying to generate localized energy 
spots be means of the FW, and not dealing with ultrasound imaging. 



transducer itself is of key importance, and perhaps the major 
issue to be considered at first in any particular application! 49l. 

As mentioned earlier, the problem of the rings definition 
implies in its turn the sampling procedure. An example of 
this process is shown in Figs. |6] and [7| both corresponding 
to the FW envelope of Fig. [8] described in Case 1 of Sec. 5. 
Here the theoretical FW profiles of amplitude and phase at the 
aperture (i.e., at z = 0) are sampled by adopting a constant- 
step approach. The sampled values are then used as the final 
amplitudes and phases of the sinusoidal excitations for each 
one of the transducer rings. The sampling process has other 
alternatives to be adopted, as for example a non-constant 
step spacing. However, we prefer to adopt in this paper the 
rather simple approach of constant sampling, leaving other 
possibilities for future work. In other words, once the values 
for d and A^^ have been assigned, we use them for a radial 
sampling of the theoretical FW patterns at z = 0: That is to 
say, each sample value corresponds here to the center of each 
ring (i.e., the center of the segment associated to the ring width 
d). 

V. Simulation OF FWs: Results 

To further demonstrate the possibilities offered by the ultra- 
sound frozen waves, in this Section we present four different 
simulations of ultrasonic FWs in a water-like medium, with 
high transverse localization [at this stage, we prefer to test 
simple arbitrary patterns]. As mentioned at the beginning. 



Fig. 9. The FW pattern, corresponding to the previous Figure, that is, to Case 
1 of Sec. 5, obtained by our simulated experiment. The parameters are set as 
in Fig.3: Namely, = 12; R = 35 mm; d = 0.3 mm; = 0.05 mm; 
Nr = 101; L = 30 mm and /o = 2.5 MHz. 



we assume in this paper a homogeneous medium, with no 
attenuation effects, operating in a Hnear regime; the sound 
speed being c = 1540 m/s. In all cases, the same annular 
aperture is used, with Nr = 101 rings endowed with width 
d = 0.3 mm and kerf of A^^ = 0.05 mm. The operating 
frequency is fixed at /o = 2.5 MHz, which is an adequate 
value for use in Medicine, corresponding to a wavelength of 
about A ^ 0.62 mm. The only parameters varied during 
our simulations, apart from the FW pattern itself, are the 
maximum allowable distance or field-depth L, and the value 
of N determining the number, 2N + 1, of Bessel beams in 
Eq.Q. 

The double Figures associated with this Section 5 (namely, 
Figs.8-9, 10-11, 12-13, 14-15) depict, first, the theoretical 
pattern to be constructed (that is, a 3D plot of the chosen FW); 
and, second, the result of the corresponding impulse response 
simulation; respectively. The sampling frequency used in the 
IR method is fs = 100 MHz. 

Let us comment about the approximate generation — by our 
simulated experi-ment — of the chosen FWs. Besides playing 
with the value of N, which enters relation ([16]), we used the 
tool of pushing the value of a little bit up, for a better 
reconstruction of the desired ideal intensity patterns |F(z)p. 
In fact, we are still using for the aperture the same size, R = 
35 mm, that is to say a diameter = 70 mm. Then, if 
is moderately increased, say for instance from 9 to 12, this 
trick will moderately enhance the reconstruction of the pattern. 
Care should be taken, however, when the maximum range L 
is augmented, because the pattern would start distorting. Such 
an effect can be appreciated in Fig [13] where the FW peak, 
near 2: ^ 55 mm, starts to loose amplitude. 

Another issue we like to comment about changing N, is 
the variation in the size of the spot-radius of the FW given 
by Eq.(|9]). However, in the present context this has not 
much concern, because all spots are below 1 mm. It is also 
interesting to compare the long depth of field of the generated 
FW, with that of a gaussian bearrj^ which in the best case 
with the current setting would be z^j^if ^ 4 mm. 

Because of the computer time employed by the IR simula- 
tions, the details in the corresponding spatial grids (p, z) had 

^The diffraction length is give by z^if = V^ko . See pp. 8-9 of Ref.fTl. 



Fig. 10. The theoretical FW pattern chosen in Case 2. Settings: = 12; 
R = 35 mm; d = 0.3 mm; = 0.05 mm; Nr = 101; L = 40 mm; 
and fo = 2.5 MHz. 




Fig. 11. The FW pattern, corresponding to the previous Figure, that is, to 
Case 2, obtained by our simulated experiment. The settings are the same as 
in the previous Figure. 



to be reduced w.r.t. the ideal patterns. Then, suitable intervals 



of A, 



0.25 mm and = 0.15 mm have been selected 



for the simulated plots. Also, due to the adoption of colors 
for the matlab graphs (absent, however, in the version printed 
on paper), an effect was added to enhance the visibility of the 
smaller details of the patterns. 



A. Case 1 

As a first choice, the FW to be reproduced consists in two step- 
functions with different amplitudes. The maximum distance 
selected for this case is L = 30 mm, the number of the Bessel 
beams to be superposed being 2N + 1 = 25 . The the spot 
size is approximately 2Api ^ 0.72 mm. The corresponding 
envelope function Fi{z) is therefore 



0.5 for li<z<l2 
1 for ^3 < z < ^4 



(18) 



with h = L/6, I2 = 2L/6, ^3 = 3L/6 and h = bL/6; the 
field being zero elsewhere. 

The theoretical pattern is shown in FigjSj while the IR 
simulation is shown in Fig|9] Notice how the peaks and the 
valleys of the ideal FW (which corresponds to = 12) are 
clearly emulated by the results of our simulated experiment 
depicted in Fig|9] Only in the region near z ^ b mm the 
simulation begins to deviate from the theoretical behavior. 







Fig. 12. The theoretical FW pattern chosen in Case 3 (see the text). Settings: 
N = 11; R = 35 mm; d = 0.3 mm; = 0.05 mm; Nr = 101; 
L = 60 mm, and /o = 2.5 MHz. 



Fig. 14. The theoretical FW pattern of Case 4. Settings: = 8; R = 
35 mm; d = 0.3 mm; = 0.05 mm; Nr = 101; L = 80 mm, 

/o = 2.5 MHz. 




Fig. 13. The FW pattern, corresponding to the previous Figure, that is, to 
Case 3, obtained by simulated experiment. The settings are as above. 



B. Case 2 

The second pattern, selected for the FW to be created, con- 
sists in a concave-shaped region, whose envelope function is 
therefore 



F2(z) = 1 + 2.5 



{z-h){z-l2) 



for h<z<l2 (19) 



with li = 2L/8 and I2 = 6L/8. The maximum distance 
adopted in this case is L = 40 mm, and the number of 
the superposed Bessel beams is again 2N + 1 = 25. This 
corresponds now to a transverse spot of 2Ap2 — 0.82 mm. 



Figures [TO] and [TT] show the corresponding theoretical and 
simulated ultrasonic field, respectively. Again, a good match 
is obtained between the theoretical and the simulated FW, 
notwithstanding the rather strange look of the pattern. 

C. Case 3 



This case corresponds to the field depicted in Figures [12] and 
T3\ Here we purpose to create two small convex, or peaked, 
field regions, followed by a constant-field finite region (a 
"step"); the field being assumed to be zero elsewhere. The 
envelope function employed this time is therefore given by 
(the field, let us repeat, being otherwise zero): 



Az-h){z-l2) 

Az-h){z-le) 
1 



for h < z < I2 
for I3 < z < I4 
for h < z < Iq 



(20) 




Fig. 15. The FW pattern, corresponding to the previous Figure, that is, to 
Case 4, obtained by our simulated experiment. The settings are as above. 



1.5L/8, /3 = 2L/8, h = 3L/8, 
7.5L/8. The obtained spot-size is 



with h = 0.5L/8, I2 

= 5.5L/8 and Iq 
2Ap3 ^ 0.98 mm. 

As we stated before, we can see the effect on the second 
peaked region produced by augmenting the number of the 
beams, which from Eq 16 should be now N' = 11 instead of 
N 7. Also, the field-depth is now larger, that is, L = 60 mm. 

D. Case 4 

The last example corresponds to an exponentially growing FW, 
with a maximum distance now even larger, of L = 80 mm. 

This is a rather important example, since it is directly related 
to the generation of a FW in an attenuating medium. The 
associated envelope function F4{z) is now 



Fiiz) 



e L 
~5~ 



for h < z <l2 



(21) 



with h = 0.5/8 and h = 7.5L/8. 

In order to enhance the fidelity in the reconstruction of the 
ideal ^4(2;) envelope, and at the same time to reduce the effects 
produced by using augmented and L values, the number 
of the superposed Bessel beams is now taken to be 17 (N = 
8). The transverse spot-size this time results to be 2Ap4 ^ 
1.36 mm. 



Notice how the ultrasonic field in Fig 15 simulates in a 
pretty good way the theoretical 3D plot of the FW in [14 
Good results of this type encourage to generate frozen waves, 
even in attenuating media. 
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VI. Some Conclusions and Prospects 

In this paper we have shown how adequate superpositions 
of zero order Bessel beams can be used for constructing 
ultrasonic wavefields with a static envelope, within which 
only the carrier wave propagates; for simplicity, we have 
here assumed a water-like medium, disregarding attenuation. 
Indeed, we have demonstrated by simulated experiments, via 
the impulse response method, that suitable sets of annular 
apertures can produce ultrasonic fields "at rest": which have 
been called ultrasonic frozen waves (FW) by us. Such FWs 
belong to the realm of the so-called localized waves (LW), 
or Non-diffracting Waves, which are soliton-like solutions to 
the ordinary linear wave equation which since long time have 
been shown theoretically and experimentally to be endowed 
with peak- velocities V ranging from to oo Q, (21. 

The FWs, among the sub-luminal (or rather sub- sonic, in 
our case) LWs|16l, are the ones associatedQ, O with F = 0. 

An important characteristic of the FWs is that they can be 
constructed, within a prefixed spatial interval < z < L, with 
any desired longitudinal intensity shape, maintaining at the 
same time a high transverse localization. This has been verified 
through the size of the intensity spots of the generated FWs, 
which have almost reached the diffraction limit A. All this 
can be monitored by appropriate selection of the parameter Q. 

We also pointed out how — when adopting for instance sets 
of annular radiators — a proper generation of ultrasonic FWs 
with a given frequency requires careful attention to the dimen- 
sioning of the annular transducers. Namely, we have discussed 
how the transducer radii (i?), the operating frequency (/o), 
and the ring dimensions (width= d and kerf= A^) affect the 
generated FW fields. We have shown, in particular, how a rise 
in frequency allows improving the fidelity of the generated 
FWs while keeping the emitter radii R unchanged, and just 
working on (augmenting) the parameter N\ even if such an 
increase of N will obviously impose sharper requirements on 
the emitter ring-elements. 

Associated to the problem of defining the ring sizes, another 
important issue we have pointed out is the required sampling 
process of the FW patterns at the aperture: Something that 
one must bear in mind when designing the annular arrays for 
the considered FW. To address this point, we explored in 
this article a simple constant-step sampling approach, leaving 
other alternatives for future work, in which we also intend 
to consider pulsed exitationsfll. |[T7 l as a possible way for 
surmounting the known side-lobe problem, met with the use of 
Bessel beams for focusing ultrasonic energy in CW (continuos 
wave) mode. 

As a merely complementary issue, we discussed incidentally 
the effect of raising even more the value of [instead of 



sticking to the value yielded by Eq. (16)] for getting in more 
detail the ideal envelope F{z) (see Section 2), neglecting for 
the moment the transverse- spot size changes. 

To spend a few more words on the practical realization of 
acoustic FWs, let us recall that significant efforts from the 
point of view of the hardware could be requested. We feel 
that the main difficulty is linked however with proper choice. 



and development, of the ultrasonic transducers themselves: 
Namely, of transducers with suitable dimensions and A/4 
adaptation layer(s), as well as with a low parasitic inter- 
element cross-coupling (e.g., lower than — 35dB). According 
to us, with one such transducer in the hands, acoustic FWs 
can be experimentally tested. Further work and attention will 
have to be devoted, in any case, to ease up the mentioned 
technological complexities, and costs, which in principle can 
be involved in the construction of the suitable CW driving 
electronics for the required annular array. Without forgetting 
that we confined ourselves, for obvious reasons, to study 
the region located between the aperture and the propagating 
medium, disregarding at this stage possible issues referring to 
the transfer functions associated with each one of the annular 
electromechanical elements: Something that could have indeed 
its role when aiming at an efficient implementation with pulsed 
HV exitations. In the present work, the latter problems do not 
show up, in practice, since we refer to operation in the CW 
regime; so that all the ring elements, plus the channels of the 
electronic front-end (LECOUER tm), are tuned at the resonant 
frequency of the array piezoelectric dye. 

Many applications of the localized FWs in various sectors 
of science are possible. In the case of ultrasonic FWs, let 
us mention for example new kinds of acoustic tweezers, bis- 
touries, and other possible medical apparatus for the treatment 
of affected tissues. With respect the last point, we believe 
that FWs can offer for example an alternative to HIFlQ 
since they allow an extended control and modeling of the 
longitudinal field envelope. A further advantage of acoustic 
FWs is that they don't have to deal with the drawback, of 
relative small treatment volume versus relatively large access 
window, usually met L48J by the HIFU techniques. 
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